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Abstract 

We compute the shear viscosity in the 0{N) model at first nontrivial order 
in the large expansion. The calculation is organized using the expansion 
of the 2PI effective action (2PI-1/A^ expansion) to next-to-leading order, which 
leads to an integral equation summing ladder and bubble diagrams. We also 
consider the weakly coupled theory for arbitrary A^, using the three- loop ex- 
pansion of the 2PI effective action. In the limit of weak coupling and vanishing 
mass, we find an approximate analytical solution of the integral equation. For 
general coupling and mass, the integral equation is solved numerically using a 
variational approach. The shear viscosity turns out to be close to the result 
obtained in the weak-coupling analysis. 
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1 Introduction 



Due to the recent progress in relativistic heavy ion colhsions and cosmology, detailed 
theoretical investigations of the dynamics of quantum fields out of equilibrium have 
become an active subject of research. A successful approach to solve the dynamics of 
quantum fields far from equilibrium as well as the subsequent stage of equilibration 
and thermalization makes use of the so-called two-particle-irreducible (2PI) effective 
action (see Refs. 0-^2] for recent nonequilibrium applications). 

The final stages of thermalization in systems out of equilibrium with conserved 
quantities can be described by hydrodynamics, characterized by an equation of state 
and transport coefficients. Recently, we investigated the connection between the 2PI 
effective action and transport coefficients for a variety of field theories and found that 
the lowest nontrivial truncation of the 2PI effective action determines correctly trans- 
port coefficients in a weak coupling or expansion at leading (logarithmic) order 
[T^.^ We emphasized that this result provides an important benchmark to validate 
commonly used truncation schemes for nonequilibrium quantum field dynamics.^ 

In applications to heavy ion dynamics, successful hydrodynamical descriptions of 
heavy ion collisions are so far based on ideal hydrodynamics, assuming infinitely fast 
thermalization and vanishing transport coefficients ^Tj. The extension to nonideal 
hydrodynamics is nontrivial, but the effects of viscous corrections are currently un- 
der investigation ^H]- In order to make further progress, it is crucial to know the 
magnitude of transport coefficients quantitatively ^H] ■ 

Transport coefficients in relativistic plasmas at high temperature can be computed 
following different approaches. The first complete calculations in hot gauge theories 
were done using kinetic theory: to leading logarithmic order in the weak coupling 
expansion to full leading order pm i21j. and in the large Nf limit [22] • Using field 
theory techniques, the shear and bulk viscosities were obtained in a single-component 
scalar field theory with cubic and quartic interactions through the summation of 
an infinite series of ladder diagrams [22] (for a more concise analysis of the shear 
viscosity, see [21])- In Ref- [2H] a simple and economical way was presented to carry 
out the sum of ladder diagrams contributing to the shear viscosity to leading order 
in a quartic scalar theory and to the shear viscosity and electrical conductivity in 
the leading-log approximation in (non)abelian gauge theories. The Ward identity in 
the calculation of the electrical conductivity in QED was studied in Ref. [2^], while 
an alternative diagrammatic approach employing a dynamical renormalization group 

^For the bulk viscosity it is necessary to go to higher-order truncations |15| . 
■^In this context it would be interesting to compute transport coefficients within the so-called 
2PPI effective action approach \W\ . 
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Figure 1: Typical skeleton diagram that contributes to the shear viscosity in the 
0{N) model at first nontrivial order in the large limit. The black dots with tiny 
lines represent the external bilinear operators. 

to study the conductivity was presented in Ref. ^7\. The prospects of extracting 
transport coefficients nonperturbatively using lattice QCD have been discussed in 
Ref. [2H] and first results have been obtained [22] (see also ISO])- In strongly coupled 
super symmetric Yang-Mills theories, the shear viscosity has been computed with the 
help of the AdS/CFT correspondence [SI]. Finally, the shear viscosity has been 
computed in the hadronic phase using phenomenological models [221 122] • 

In this paper we consider the 0{N) model in its symmetric phase and compute 
the shear viscosity to first nontrivial order in the expansion. The 0{N) model is 
widely used when applying methods and techniques in thermal and nonequilibrium 
field theory, because without running into the issue of gauge invariance, it still leaves 
a nontrivial problem to solve. Moreover, it acts as a low energy effective description 
of QCD (for = 4) and is frequently used in early cosmology (inflation, reheating). 
The large expansion offers the possibility to explore the behavior of transport co- 
efficients outside the weak-coupling domain and in combination with the 2PI effective 
action, it naturally leads to the inclusion of the required medium effects in the single 
particle propagators. We use a field theory approach to compute the shear viscosity. 
As we will show below, a typical ladder diagram that contributes at leading order 
is shown in Fig. ^ We note that the lines in this diagram are dressed propagators, 
so it should be regarded as a skeleton diagram. As we have recently shown the 
summation of these classes of diagrams is neatly organized using the 2PI effective 
action. In particular, the 2PI-l/iV expansion to next-to-leading order (NLO) [3J [H] 
sums precisely the relevant set of diagrams. 

The paper is organized as follows. In Section |2] we describe the 2PI effective action 
for the 0{N) model and write down the ensuing integral equation at next-to- leading 
order in the expansion. We argue that it sums all necessary diagrams that 
contribute to the shear viscosity at leading order. In the next Section we study the 
various elements that appear in the integral equation and which are essential for its 
calculation: the single particle spectral function, the gap equation for the mass and its 
renormalization, the auxiliary correlator summing the chain of bubbles diagrams, and 
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the thermal width. In Section^ we derive a compact expression for the shear viscosity 
in terms of an effective vertex summing the ladder diagrams. Section El is devoted to 
the integral equation for the effective vertex. We show how to cast the problem of 
solving it into a variational one, better suited for numerical analysis. In Section 6 we 
study the weak-coupling limit of the integral equation and find an analytical result in 
the limit of ultrahard momentum. In Section [7| we solve numerically the variational 
problem and present the results for the shear viscosity. The final Section is devoted 
to the conclusions. In Appendix A we consider the weakly coupled 0{N) model for 
arbitrary using the three-loop expansion of the 2PI effective action. 



2 2PI-1/A^ expansion 

We consider a real scalar A^-component quantum field (pa = 1, • • • , A^) with a 
classical 0(A^)-invariant action, 



The mass parameter mg and coupling constant Aq are bare parameters. To the order 
we are working, field renormalization is not necessary. We use the notation 



dx^ / (fx. 



(2) 



where C is a contour in the complex-time plane. We will work in the imaginary time 
formalism, so below we specialize to the Matsubara contour, running from to —i/T. 

According to the Kubo formula, the shear viscosity can be obtained from the slope 
of a spectral function at zero frequency. 



LA. 



?o=o 



where 



(3) 
(4) 



with TTj,- the traceless part of the spatial energy-momentum tensor. 



TTijix) = di(j)a{x)dj(j)aix) - -6ijdk(j)a{x)dk(j)aix) . 



(5) 



It is well known from weak coupling studies that a one loop calculation of the 
shear viscosity is incorrect, since diagrams that appear at higher order in the loop 
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expansion contribute at leading order and the computation must be carried out using 
dressed propagators [231 For a single- component scalar field (A^ = 1), Jeon 

showed that these higher order diagrams are ladder diagrams, where the rung in the 
ladder is a single bubble [22] • Due to the kinematical configuration, the propagators 
on the side rails suffer from pinching poles. As a result the thermal width, which 
is determined by the imaginary part of the self energy, has to be included in these 
propagators, while the real part is subleading and can be neglected (at weak coupling). 
The picture is quite similar for gauge theories at leading logarithmic order in the 
weak coupling limit [23 120] • The ladder series is conveniently summed through an 
integral equation for an effective vertex [221 12S Ei] • We have shown recently that 
this integral equation in the required kinematic configuration appears naturally in the 
2PI effective action formalism ^3]. In the expansion, in which we are interested 
here, the presence of pinching poles also leads to higher order diagrams in the loop 
expansion contributing at leading order, and the 2PI effective action formalism neatly 
organizes this calculation as well, as we describe now. 

The 2PI effective action is a functional of the time-ordered two-point function 
Gab{x,y) = {Tc(f)a{x)(t>b{y)) and can be parametrized as [211123 

T[G] = ^TrlnG-i + ^TtG,\G - Go) + T^IG]. (6) 

Throughout we consider = 0. The classical inverse propagator iGg ^ is given 

by 

^Go^abi^, y) = - + ml] SabScix - y), (7) 

with Sc{x — y) = Scix'^ — y°)5(x — y). T2[G] is the sum of all 2-particle irreducible 
(2PI) diagrams with no external legs and exact propagators on the internal lines. 
Extremizing this effective action, 

^rl'^l -0, (8) 



6Gab{x,y) 

leads to a Dyson equation for the two-point function 

G~Jix,y) = G,^,,ix,y)-J:a,ix,y) (9) 

where the self energy 

^ab{x,y)^2z j!^f\ (10) 
dGab{x,y) 

depends on the full propagator G. 
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Figure 2: Integral equation for the 4-point function. 



Within the 2PI formahsm, a 4-point vertex function T^'^^ appears, defined as the 
usual connected 4-point function with the external legs removed. This 4-point func- 
tion obeys an integral equation that can be obtained using standard functional rela- 
tions jMl (see Fig. E)) 



^ab;cdi^^ 3;', y') = Aab;cd{^, V, x', y') 

Aab;ef{x, y] W, z)Gee'{w, w')Gff\z, z')T''^y^^^{w' , z' ] x' , y'). (11) 



^(4) 



The kernel or rung A is defined as 



Kb;cd{x,y;x\y') 



, 5'Lai,{x,y) 
'5G,d{x\y'y 



(12) 



Semicolons separate indices with a different origin |14j . 

The relations so far are general. We now specialize to the 2PI-1/A^ expansion to 
NLO. The 2PI part of the effective action can be written as |31in] 



NLOr 



where (see Fig. EJ 



T\^[G] 



4:\N 



Gaa{x,x)Gbb{x,x), 



pNLO 
J- 2 



[G]=-TrlnB, 



(13) 

(14) 
(15) 




Figure 3: Contributions to the 2PI effective action in the 0{N) model in the 2PI-1/A^ 
expansion at LO and NLO. Only the first few diagrams at NLO are shown. 
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Figure 4: Contributions to the self energy at LO and NLO. The auxihary correlator 
D is represented by the dashed line. 



>-< = X^X3-< 



Figure 5: Integral equation for the auxiliary correlator D. 



with 



B(x, y) = 5c{x-y)- ^n(x, y). 



The function B depends on the single bubble diagram defined by 

n(a;,2/) = -]^Gab{x,y)Gab{x,y). 

Expanding the logarithm in Eq. (fT3j) generates the closed chain of bubble diagrams 
as in Fig. El The self energy follows from Eq. We write S = E^^ + 2^^°, with 



(16) 



(17) 



i^5abGcc{,x, x)5c{x - y), 



-ah y) 

^^^°{x,y)=-Gab{x,y)D{x,y), 
shown in Fig. HJ Here we introduced the auxiliary correlator 

D{x,y) = '^B-\x,y). 

From the identity B~^B = 1 it follows that the auxiliary correlator obeys 

Scix-y) + Ju{x,z)D{z,y) 



(18) 
(19) 

(20) 
(21) 



which is depicted in Fig. For the rung, we write A = A^*-* + A'^^'^'-', and find (see 
Fig. ED 



Aa6°cd(^. 2/; 2;', y') 



iXo_ 
'3N 



5ab5cd5cix - y)5c{,x - y')5c{x' - y), 



Kb;cd{x^ 2/; a;', y') = - [SaJbd + SadSbd D{x, y)Sc{x - x')6c{y - y') 

+2Gab{x, y)D{x, x')D{y, y')Gcd{x', y'). 



(22) 



(23) 
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X ^ J _ 

a) b) c) 

Figure 6: 4-point kernel in the 2PI-l/iV expansion of the 0{N) model at LO and 
NLO. 



Finally, we specialize to the 0(A^)-symmetric case, take Gab{x,y) = 6abG{x,y) 
and T,ab{x, y) = 6ab^{x, y), choose the Matsubara contour and write the equations in 
momentum space. The final set of equations then reads 



with 



u;'^ + p^ + ml + j:{py ^^^^ 



S"°(P) = y^G(/0, (25) 
SNLO(p) = G{P + K)D{K), (26) 



and 



^<^>= -3iV/Ao + n(P) - 

Tl{P) = -j^^G(P + K)G(K)^ (28) 
The 4-point function obeys 

rit,(i?, K) = Kb-AR. ^) + ^ ^ ^«'';^/(^' P)G\P)T%,{P, K), (29) 



with the kernel 



A^;,%(i?,P) = -^5a.5c., (30) 



^ai?d{R, P) = %chd + 5adhc] D{R - P) 

+28ab5cd ^ G{R - L)D\L)G{L - P). 



(31) 



Here we used the notation P = {iuOn^ p), where Un = 27rnT (n G Z) is the Matsubara 
frequency, and 



In the remainder of this section, we argue that the set of Eqs. (l24|) - (|3T|l sums all the 
diagrams contributing to the shear viscosity at leading order in the expansion. 
First we note that the one loop contribution to the viscosity is proportional to A^^. 
Here, one factor of arises from the group indices running in the loop, the other 
originates from the pair of pinching poles in the loop. Pinching poles are screened 
by the imaginary part of the retarded self energy E, which appears first at NLO (the 
leading order contribution to S is real). Since pinching poles are sensitive to the 
inverse of the imaginary part, each pair of pinching poles gives a contribution that 
scales as A^. Therefore the shear viscosity takes the form^ 



where m/j is the renormalized mass in vacuum and A(/i) is the renormalized coupling 
constant at the scale n. The shear viscosity is independent of /i. It is now straight- 
forward to identify which diagrams contribute to the shear viscosity at leading order. 
The auxiliary correlator D is proportional to see Eqs. (j2Zll2HI)- Starting from 
the naive one-loop expression, adding a vertical D correlator as in Fig. yields 
one extra pair of pinching poles that cancels the explicit 1/N from the D correlator. 
Therefore all vertical line insertions contribute at the same order. An insertion of the 
box rung (see Fig. IHt) results in two additional D correlators, one additional pair of 
pinching poles, and one additional closed loop over the group indices. Therefore also 
all box rung insertions contribute at the same order. Finally, the lowest order rung 
(see Fig. IHK), which would generate a string of bubbles, does not contribute to the 
shear viscosity. This was shown by Jeon for the weakly coupled single compo- 
nent case and will be confirmed below. Rungs that are down by 1/A^ have additional 
D correlators but no compensating additional pinching poles or closed loops. Some 
examples of subleading rungs are shown in Fig. [7| We find therefore that the typical 
diagram contributing at leading order is as shown in Fig. ^ Throughout the paper 
we neglect subleading powers of A^ as indicated in Eq. (j221)- 

Before analyzing the integral equation, which sums the required contributions, 
several elements are needed in detail. We study those in the next section. 

■^Note that the ratio of the shear viscosity and the entropy density s is proportional to N, and 
therefore always far above the lower bound conjectured recently, 77/s > l/An |36j . 





(32) 




(33) 
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Figure 7: Examples of rungs that may contribute at next-to-leading order to the shear 
viscosity in the expansion. 



3 Quasiparticles 

In this section we study the single particle propagators that are dressed to account 
for plasma effects, which are crucial to the calculation of transport coefficients. A 
convenient way to study these effects is through the use of the single particle spectral 
density. 

Analytically continuing the euclidean propagator (|24j) to real frequencies yields 
the retarded and advanced propagators, 



Gr{P) = G{iujn ^ / + zO, p) = G*AP) 

1 



-pI + p2 + + ReS«(P) + ^lmS^(P) ' 
from which the spectral function is obtained as 

p{P) = -I [Gr{P) - Ga{P)] 

-2ImS^(P) 



(34) 



[pI - p2 - - ReT.R{P)f + [ImS^(P)]'' 



(35) 



The leading part of the self energy is of the same order as the mass in the 
expansion and hence it cannot be neglected. The real part of on the other hand 
can be dropped since it is suppressed by compared to the leading piece. ^ Since 
ImS/j is proportional to the spectral density generically has a simple on-shell 
form in the large N limit, 

p(P) = 27rsgn(/)5(p2 _ + 0(l/iV), (36) 

with 



cUp = Vp' + M\ M^ = ml + S^o. (37) 

However, when pairs of propagators with pinching poles are present, the imaginary 
part cannot be neglected as it separates the poles in the complex-energy plane that 

■*This is the reason field renormalization does not enter in this problem. 
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approach the real axis from below and above [231 123- In this case, the full propagator 
G{P) has a cut along the whole real axis since ImE/j(P) is nonvanishing for all 
real p^. Moreover, the product of retarded and advanced propagators is directly 
proportional to its inverse, 

Gn{P)GAP) = _ _ ^2 _ ReS^(P)]^ + [ImS«(P)]^ = -2Ims!(P) ' ^^^^ 
In the limit of large we may use Eq. and this expression reduces to 

Gn{P)GA{P) = ^^+0{l), (39) 

with Fp the (on-shell) thermal width defined from the imaginary part of the retarded 
self energy as 

ImS^(P) 



r 



P pO 



(40) 



In the remainder of this section we study the corrections introduced by the self energy. 
The real part of S leads to a gap equation for the mass which requires renormalization 
of both the mass parameter and coupling constant. The imaginary part is computed 
in terms of the auxiliary correlator D, which we work out in detail. 



3.1 Gap equation and renormalization 

Since the real part of S^'"^ can be systematically neglected, the gap equation for the 
mass parameter M reads 

M' = ml + ^^^G{P), (41) 

with 

- (42) 

This gap equation is divergent and in order to renormalize it^ we also need the integral 
equation for the 4-point function, Eq. ()29p. at lowest order. In this approximation 
the 4-point function is momentum independent and Eq. can be solved as 

k) = 5..5..r(^), -L = + n(o), (43) 



'See Refs. (33 |2H] for recent studies of renormalization in the 2PI effective action formalism. 
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with the single bubble at zero momentum 

n(o) = -y^G^(P). 



(44) 



We note that the equation for the 4-point function is identical to the equation for the 
auxiliary correlator D{P) at zero momentum. Therefore the renormalization carried 
out to obtain a finite gap equation renders D{P) finite too. 

In order to regulate the divergent integrals we use dimensional regularization in 
3 — 2e dimensions. Both Aq and F^^^ have dimension jj?^. We introduce a dimensionless 
coupling \r via 

^^^'^ (45) 



r(4) 



3A^ 



such that the equations to renormalize read 



rrin 



(46) 
(47) 



Since renormalization can be carried out at zero temperature, we only compute the 
bubble in the vacuum 



1 



2e 



Ao 



1 



967r2 



- + ln(47r) 



7ij + 21n-^ 



(4^ 



where mR is the renormalized mass at zero temperature. Introducing the running 
coupling constant in the MS scheme A(yu) via 



A(/i) 



we find 



Ao ^ 967r2 



1 



+ ln(47r) - 7£; 



In 



/i 



A(/i) 487r2 niR 



(49) 



(50) 



The running coupling A(/i) obeys the usual renormalization group (RG) equation with 
the correct (5 function for the 0{N) model in the large limit, /3(A) = A^/(487r^). 
Returning to the gap equation (jTTj) at zero temperature, we find 



m 



R 



rrin 



-2e 



m 



R 



6 167r2 



+ ln(47r) - + 1 + 2 In 



rriR 



(51) 
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Figure 8: Finite temperature mass M/T as a function of the renormalized mass at 
zero temperature rrtji/T for 3 different values of A(T). 



Combining Eqs. (|471 EHl we arrive at the renormahzed gap equation at finite 
temperature 



mi 



967r2 



m 



R) 



1 + ln 



m 



R 



m. 



^^^^^ 



6 



niujr 



UJr 



(52) 



where n{uj) = l/(e^^'^ — 1) is the Bose distribution function. The solution of this gap 
equation is independent of the renormalization scale /i. A numerical solution of the 
gap equation for three values of A(/i = T) is shown in Fig. |H1 

There is one additional issue related to renormalization that needs to be consid- 
ered. This scalar theory has a Landau pole at the scale A^, = fxe^^'^ ^ where the 
running coupling constant diverges. One has to require that all physical energy scales 
involved in the problem, i.e. tur, M, and T, are much smaller that the one associated 
with the Landau pole. This imposes a restriction on which values of A(/i) and m^j can 
be considered. Taking, for example, T to be the largest scale and demanding that 
A^/T > 40, we find that A(/i = T) < 487rVln40 ^ 128. 



13 



3.2 Auxiliary correlator 

Before proceeding with a calculation of the thermal width, we need to consider the 
chain of bubbles summed by the auxiliary correlator D. The single bubble is given 
by 

n(P) = -f ^ G(P + K)G{K), (53) 
which after performing the Matsubara sum reads 



1 



+ [1 + n(cjk) + n{ur)] ■ ■ \ , (54) 

where r = k + p. In order to separate the logarithmic divergence in the zero- 
temperature contribution we write 

n(P) = no(P) + nT(P), no(P) = no(o) + n'o(P), (55) 

where no(P) is the part without distribution functions. The divergent contribution 
no(0) is similar to the one computed for the renormalization of the coupling constant, 



no(0) = Q + ln4vr - + 2 In ^ ) , (56) 



while the remainder is finite and can be evaluated as 



with P^ = cj^ + p^. The easiest way to arrive at the above result is to go back to the 
original four- dimensional euclidean integral. 

We are interested in the retarded bubble, obtained by analytical continuation 
iujn — ^ p° + iO. For the finite contribution to the real part we write 

Ren^(P) = Reno^(P) + Renf (P), (58) 

with 

-m 



Ren;«(P) = "Y^jl + \Hs- am') + &{-s)] (3{P) In 



1 + /5(P) 



-e(4M2 - s)e(s)P(P)arctan-^^|, (59) 
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where s = pg — p^, and 

m = 

and 

Ren^(P) = 



4M2 

2 9 ' 



B{P) 



- 1, 



N 



°° k 
dk — n(ix'k) In 



{k + p+){k + p.) 



{k - P+){k - P-) 



(60) 



(61) 



with p± = I [p ± //3(P)]. The remaining integral can be done numerically. 

The imaginary part can be written in terms of single particle spectral functions 

as 

lmUn{P) -- 



N 

T 



K 



where 



p{P + K)p{K) [n{k^) - n{p^ + k^)] 
d^k 



K 



(62) 

(63) 



and can be evaluated completely 



Imnij(P) = -0(s - 



N 
'32^ 



2T 1 — p-P+l'^' 

/3(P) + -ln \ ,^ 
p 1 — e P-'^ 



NT 1 - e~P+/'^ 
.e(-s)^ln^^— ^. 



(64) 



withpi = l[p°±pf3{P)]. 

The chain of bubbles is summed by the auxiliary correlator D in Eq. ()27|) . The 
divergent piece of the single bubble no(0) is absorbed by coupling constant renormal- 
ization. The renormalized expression for reads 



-3N 



In 



_A(/i) 967r2 _ 



U',{P)+Ut{P). 



(65) 



The auxiliary correlator is renormalization group invariant. Retarded auxiliary prop- 
agators are obtained as Dr{P) = D{iujn p^ + iO, p) = D\{P). The spectral density 
is 



Pd{P) = -i [Dr{P)~Da{P)] 



-2Imnfi(P) 



3A^/A(M) - Reni(P) + [\mI{R{P)X 



(66) 
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Figure 9: Spectral density Np£,{uj, p)/A^(M) as a function of u/T for fixed p/T = 2, 
mji/T = 1 and three values of the coupling A(/i = T). The vertical line indicates the 
lightcone. 



in terms of the RG invariant coupling 



1 1 1 

In^. (67) 



A(M) A(/i) 967r2 jj,^ 

The spectral function is shown in Fig. IHl for a typical choice of parameters. We 
scaled out the trivial N/\^{M) dependence. Imn/j(cij,p) and therefore p£)(u;,p) are 
nonzero below the lightcone (u;^ < p^) and above threshold (cj^ > + 4M'^). We 
note that for larger coupling constant the contribution below the lightcone diminishes. 
Fixing mji and increasing A(T) results in a larger value for M (see Fig.|S)). As a result 
the contribution above threshold starts at larger u when increasing A(T). 

3.3 Thermal width 

The thermal width is given by 

^ ^ ImSflXP) ^ ImEnjup, p) ^^^^ 
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1.5 



p/T=5 — 



\ 

in/r=0 
m/r=0.5 



.-. m/r=i 

_. m/r=2 




20 40 60 80 100 120 

Figure 10: Width Fp/Foo as a function of A(T) for fixed p/T = 5 and four values of 
rriR/T. 

The LO part of S does not contribute. We write the self energy at NLO as 

SNLO(p) = _^G'(^)^(^_P)_ (69) 

After performing the Matsubara sum and the analytical continuation, we find ImS^ 
in terms of spectral functions as 



ImE^(P) = ]- f p{R)pd{R - P) [n(r°) - n(r° - / 
^ Jr 



pr '^pr j 1 



A convenient way to proceed is to introduce k = |r — p| as 

i-r+p ^ 
'\r-p\ 

where cos ^p,. = P ■ f is the cosine of the angle between p and r and 



poo I-r+p ^ 

1= / dk 5{k — \r — p\) = / dk — 5{cos9. 
Jo J\r-p\ rp 



(70) 



(71) 



2rp 



(72) 
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We perform the integral using p{R) and the 6pr integral using the delta function 
introduced above. The final result for the width then reads 

1 f°° r f+P ( 
Tp = — / dr — / dkk\ poi^^r + ^^p, k) [n{uJr) - n{uJr + uJp)] 

IGTT^pUJpJo J\r-p\ ^ 

-pni'^v - t^p, k) [n{uj,,) - n{uj,, - Up)] \. (73) 



The remaining two integrals can be performed numerically. Note that to obtain pD, 
a numerical evaluation of the real part of the bubble is required as well. 

In Fig. ^1 we show the width as a function of the coupling constant for fixed 
p/T = 5 and several choices of m^/T. In order to remove the trivial parameter 
dependence, we rescaled Fp with 

r ^mrn (74) 

the thermal width at ultrahard momentum p ^ T in the weakly coupled, massless 
limit at leading order in the expansion (see Sec. We observe that for small 
mass the dependence on the coupling constant is substantial, whereas for larger mass 
it becomes negligible. 



4 Shear viscosity 

Now we have all the ingredients necessary to compute the shear viscosity. In order to 
do the Matsubara frequency sums, it is convenient to transform the integral equation 
for the 4-point function into an equivalent one for a 3-point vertex, defined in Fig. ^2 
The corresponding integral equation is depicted in Fig. ^| We denote the effective 
vertex as Fjj ^^(P + Q, P) where P is the (euclidean) momentum on the siderails and 
Q = {iujq, 0) the momentum that enters from the left. The integral equation then 
reads 

TijAP+Q, P) = ^,abiP)+l ^ G(P+g)F,,,,rf(P+g, R)G{R)A,,,ab{R, P; Q), (75) 




Figure 11: Relation between 3 and 4-point vertex function. 
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Figure 12: Integral equation for the full 3-point vertex. 




Figure 13: The correlator (iiiji^ij) in terms of the full vertex. 

where the kernel Acd-ab{R, P] Q) for vanishing Q was presented in Eqs. fjHIH EIJ and 
the bare coupling between the scalar field and the Tiij operator follows from Eq. (jnj: 



PiPj - ^S^jp'^ 



(76) 



The correlator we need to obtain the shear viscosity takes the form of an effective 
one-loop diagram (see Fig. [T^. 



(77) 



The effective vertex can be taken of the form^ 



r.„a6(P + Q,P) = 25, 



ab 



PiPj - g'^ijP^ 



r{p + Q,p), 



which yields 



4 



= ^^p'G{P + Q)T{P + Q, P)G{P). 



(79) 



In order to do the sums over the Matsubara frequencies, we follow the technique 
presented in Ref. [22], which makes use of the following relation 



Tj]/(^a;„) = E 



cuts 



2m 



-n 



(ODisc/- J]n(^,)Res(/,. 

poles 



(80) 



^In principle, there could be an additional term proportional to Sij. However this would not 
contribute when contracted with the traceless bare vertex. 
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which requires knowledge of the analytical structure of the function to be summed. 
Because of the inclusion of the NLO contribution to the self energy in the propagators, 
the propagator G{p^, p) has a cut along the entire real axis, i.e. when \m.{jp) = 
in the complex p° plane. It follows from the integral equation Eq. ()75|) that the 
effective vertex has the same analytic structure, namely, r(P + Q, P) has cuts when 
Im(p'' + = and Im(p'') = 0. With this information we can do the Matsubara 
frequency sum in Eq. (j79|) . make the analytical continuation iUg + iO afterwards 
to arrive at the retarded function and take the limit — > (see Refs. [23 12^] for 
further illustration). Only the product of retarded and advanced propagators suffers 
from pinching poles and dominates at leading order in the expansion. Therefore 
only one particular analytical continuation of the full vertex is needed. Defining 

P(/,p)= lim Rer(p° + g° + «0,p°-«0;p), (81) 

the result for the spectral density reads 

hm p..(g°,0) = -^Nq' I pV(/)Gr(P)G^(P)P(/, p). (82) 
go^o 6 Jp 

Using Eq. (j39|) for the product of Gr and Ga and definition Q for the shear viscosity, 
we get 

--5X»''^v|^-<^°>')-S/;»v4^. (S3) 

where we used that n'{—ujp) = n'{ujp) and we defined 

V{p) = V{±ujp,p). (84) 

To arrive at our final expression for the viscosity, we proceed as in the case of gauge 
theories [23 120] and define the dimensionless quantity 

59^ x>(p) ^ ^ 

X{p) = (85) 

The shear viscosity then reads 

N f p^ N f°° p'^ 

^ = -77 / — ^'(^p)x(p) = -^TT^ / dp—n'{u;p)xip). (§6) 
15 Jp ujp SOtt^ Jq ujp 

Since the thermal width is inversely proportional to N, it follows that x{p) ^ N in 
the large N limit. 
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5 Integral equation 



As we have shown in the previous section, in order to obtain the shear viscosity, we 
need a particular analytic continuation of the effective vertex in the limit g° — and 
p° = itcUp. In this section we first explicitly write down the integral equation in this 
kinematical limit. We then show how to cast it in a form suitable for a variational 
treatment. We also point out the relation between the integral equation and kinetic 
theory. 

With momenta flowing as illustrated in Fig. El the integral equation reads 
r.,>6(P+g, P) = vlM+l ^ G{R+Q)T,^,UR+Q, R)G{R)A,,,ab{R, P; Q), (87) 
where the rung is 



Acd,ab{R, P', Q) 



An 



3iV 



Sabred + {^achd + ^adhc)D{R - P) 



-^25ab5cd X D{L)D{L + Q)G{R - L)G{L - P). 



(88) 



To arrive at a scalar equation we contract with ^^^^^^(p) and divide by a common 
factor to find 

r(p+g, p) = 1 + -1- ^ ^P2(p- r)G(p+g)r(p+g, p)G(P)A,e,,,(p, p- q). (89) 

2N TrV 

Here, the second Legendre polynomial P2{x) = {3x^ ~ l)/2 arises from the contraction 



^^",,a6(p)r.j,cd(P + g, P) = :sP r'P2{p ■ ^)5abScdT{R + g, P) 



(90) 



r+q 



p+q 




r+q 



l+q p+q 



q 




\p-l 



P I P 

Figure 14: Choice of momentum labeling. We also denote k = |r — p|. 
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We now perform the Matsubara frequency sums following again the technique in 
Ref. j2Sl- The lowest order term in the rung does not contribute since it is a constant 
and it has no discontinuity. For doing the sum for the second term in the rung, we 
need to know that the auxiliary correlator D{R — P) has a branch cut below the 
lightcone and above threshold (see Sec. 4.2). However, when using Eq. (jHUj) . we may 
assume that the cut runs along the entire real axis, since the contribution from where 
the cut is absent cancels automatically. Finally, for the last piece of the rung there 
are two sums. These can be performed sequentially, provided that the analytical 
continuation is made at the very end. The analytical structure is as follows. The G's 
that depend on R + Q and R as well as the vertex have branch cuts along the entire 
real energy axis, for the D's we may take branch cuts along the entire real energy 
axis, and the G's that depend on R — L and P — L contribute just simple poles (since 
they do not require the inclusion of the thermal width). When performing the sum 
over ui, one picks up contributions from the poles of the G's and from the cuts of the 
-D's. We found that after making the analytical continuation only the poles from the 
G"s contribute. 

From Eq. (j8ip it follows that we need the following analytic continuation 

iup — > p° — iO, iujq ^ g° + zO, iujp + iujq — > p° + g° + zO. (91) 

Preserving again only the dominant contribution with pinching poles for vanishing 
g", we arrive at 

V{P) = 1+1 ^P2(p ■ i) [ri(rO - /) - n(rO)] Gr{R)V{R)Ga{R)KR, P), (92) 
Jr P 

with 

AiR,P)=PDiR-P)+N j [n{p^ -f)-n{r^ -f)] p{P - L) p{R - L) \D r{L)\\ 

(93) 

We note that the kernel obeys A(i?, P) = -A(P, R), A{-R, P) = -A(i?, -P). 

We now use Eq. (jHHI) for the product of Gji and Ga, Eqs. dHU IS^ to introduce 
x{p), and take p° = ±ujp (both choices yield the same equation) to write the integral 
equation in the final form 

^prpX(p) =P' + l [ - ^p) - ri{r')] ^P2(p ■ f)x(r)p(P)A(P, P). (94) 

Jr ^ 

Solving for x{p)i one obtains the shear viscosity from Eq. (j86j) . 
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5.1 Variational approach 

Since the integral equation cannot in general be solved analytically, one has to rely 
on numerical methods. A convenient approach to follow is the one employed in 
Refs. |191 1^ 122] ! where the problem of obtaining a transport coefficient in kinetic 
theory from an integral equation is formulated as a variational problem. Here we 
show how to cast the integral equation in a form that is suitable for a variational 
treatment. 

After multiplying Eq. ()94|) with 

-n'M, (95) 

Up 

the integral equation can be written rather compactly as 

POO 

^{p)x{p) = S{p) + I drn{p,r)x{r), (96) 



with ^ 

^(p) = p'n'iup)Tp, S{p) = ^n'iup), (97) 

Up 

and a symmetric kernel 

nip,r) = nir,p), (98) 

whose explicit form is presented below. Since 7i is symmetric, Eq. can be derived 
by extremizing the functional 



Q[x] = dp 



1 1 f°° 



(99) 



From Eq. ()86|1 we see that the actual value of this functional at the extremum, 

Q[x = Xo.t] = \j^ dpS{p)x{p), (100) 



is directly proportional to the viscosity 



N 

^ = -T^QlX = ^ext]. (101) 



In the rest of this section we explicitly evaluate Ti.{p, r) 
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We treat separately the diagram with the single line and with the box diagram 
and write Ti = Hnnc + '^box- We start with the diagram containing the single line. 
The integral to evaluate reads 

[ [n(rO - Up) - n(rO)] "^P.ify ■ r)x{r)p{R)pD{R - P). (102) 
Jr 

We proceed in exactly the same way as in the case of the width in Sec. 13.31 and 
introduce /c = |r — p| via Eq. (ffT|l . We perform the r° integral using p{R) and 
the integral over the angle between p and r using the delta function in Eq. 1)7111 . 
Multiplying the result with Eq. (jHSj), we can immediately read the contribution to 
7i(p, r) introduced above, and we find 



'H\inc{p,r) = — — ^ / dkk P2{Zpr)\pD{^^r " ^^p, k) [n{ur - Up) - n{ur)] 

IbTT^ Up Ur J\p_r\ ^ 

-pDiuJr + Up,k)[n{ur + Up) -n{ur)] \, (103) 



where Zpr is given in Eq. ()72|1 . Using the relations 

, N / N Tn'iur) , N / N Tn'iur) 

n{Ur + Up) — n[Ur) = -, r -, r, n[Ur — Up) — n[Ur) 



1 + n{ur) + n{up) n{ur) — n{up) 

(104) 

and the fact that poi^r — ^p, k) is odd under interchange of p and r, it is straight- 
forward to verify that the result is symmetric in p and r. 

For the diagram containing the box rung, we have to evaluate 

[ [n(r° - Up) - n(r°)] [n{up - /°) - n{r^ - f)] 
Jr,l 

xP.iP ■ r)x{r)'^p{R)p{P - L)p{R - L) \Dn{L)f . (105) 

There are three angular integrations that are nontrivial. We denote the cosine of the 
angle between p and 1 as cosOpi, between r and 1 as cosOri, and the azimuthal angle 
between the p, 1 plane and the r, 1 plane as </>. The 8-dimensional integral can then 
be written as 

r\ /"OO roo roo roo /»! rl /•27r 

— — / drr^ dr^ / dlP df / dcosOpi / c/cos^w / dcf). (106) 

(27r) Jo J-oo Jo J-oo J-1 J-i Jo 

The integration over cos 6pi will be performed using the delta functions in p{P — L), 
the one over cos 6ri using p{R — L) and that over r° using p{R). The product of the 
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three spectral functions yields a set of constraints, since 

p{R)p{P - L)p{R -L) ~ y (5(r° + sicUr)<5K - /° + S2c^p_i)5(r° - /° - s^uo^.i) 

~ y (5(a;p + SiUJr + S2Ci;p_i + Ss^r-i). (107) 

Si=± 

Out of the eight combinations, only three can contribute for kinematical reasons, 
namely those corresponding to 2 ^ 2 processes. We treat these three cases separately 
and write 

Hbox = H« + + (108) 

1. (si, S2, S3) = (— , +, — ). The cosines are cos9pi = Zpi, cos9ri = z~i, where 
The constraints from the spectral functions can be satisfied provided 



/° > + 4M2, |/_| <p,r < |/+|, (110) 
where here and below 



1 / 4A/f2 

4 = -[/±mL)], /^(^) = yi-^T37^- ^^^^^ 

The only place where the angle (p appears is in p ■ f , which when expressed in 
terms of the angles we use, reads 

p • r = sin 6pi sin O^i cos (p + cos 6pi cos 6^1. (112) 

It is then also straightforward to perform the 0-integral and we find 

/ rf0P2(p ■ f) = 27rP2(cos^p0^2(cos^w). (113) 
Jo 

Multiplying the resulting expression with Eq. (|95p we can read off the first 
contribution to T-C{p,r) from the box diagram: 

/•oo /»oo 

X dl dfP2{zpi)P2{z;^) \DniL)f [n{up - /°) - n(^r - l' 

Jo JVW+iM^ 

xe(p- |/_|)e(|/+| -p)e(r- |L|)0(|/+| -r). (114) 
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2. (si,S2,S3) = (—,—,+). The cosines are cos^p; = Zpi, cos9ri = z^^ and the 
constraints are 

^g</^ P>|/+|, r>|/+|. (115) 
Therefore the second contribution reads 
^2). x_ N p r 



T^hLiP^r) = — — n'{ujp) [n{uJr - ujp) - n{ujr)] 

iZoTV^ Up UJr 

f oo 



dl df P2{zpi)P2{z-) \Dr{L)\^ [n{ujp - f) - n{uo, - f)] 
xe(p-|/H.|)e(r-|/+|). (116) 

3. (si,S2,S3) = (+,—,—). In this case the cosines are cos9pi = Zpi, cos9ri = z^i, 
with the constraints 

/^</^ P>\l+\, r>|/_|. (117) 
The third contribution then reads 
.(3) _ J^P^^n 



T^hLip^r) = '^''^^p) ['^'^^r + ^p) ~ '^(^r)] 

f oo pi 



XI dl j df P2{zpi)P2{z+) \Dn{L)f [n(/° - Up) - n{f + uj, 
xe(p-|/+|)e(r-|/„|). (118) 



Using relations ()1()4|1 and making the substitution /° in the third contribution, 

one finds that Hhoxip,^) = 'Hhox.ij.p). In conclusion, 7Y(p, r) = H{r,p), which allows 
to obtain the integral equation from the functional Q. 

5.2 Kinetic theory 

Here we briefiy mention the relation between our results and the corresponding kinetic 
theory, by analyzing the kernel in the integral equation. 

The kernel A{R,P) can be written in a form that allows for a direct comparison 
with kinetic theory. We start by changing variables L ^ P — L and find 

A{R,P)=Pd{R-P)+N J [n(/°)-n(r° -/ + /")] p{L)p{R-P+L) \Dr{P - L)\\ 

(119) 



26 




Figure 15: Scattering processes in kinetic theory, time runs horizontally. 



If we use for the spectral density the following identity 

Pd{R -P) = -2Imn^(i? - P) \Dr{R -P)\\ (120) 
where (see Eq. (jHSI)) 

ImUniR -P) = -^ [ [n{f) - n(r° - / + f)] p{L)p{R -P + L), (121) 

introduce the variable L' according to 

/ (27r)V(i?-P + L-L'), (122) 
Jl' 

and use the delta function to interchange momentum labels, the kernel can be written 

as 

A{R, P) = ^ [ [n{f) - n{n] p{L)p{L'){2'Kf5\R - P + L ~ L') 

^ J L,L' 

X \DR{R-P)\^+\DR{R-L')\^ + \DR{R + L)f^. (123) 

The second line is precisely the sum of the squares of the matrix elements for 
scattering (i?, L) — > (P, L') in kinetic theory, where the auxiliary correlator carries 
the interaction (see Fig. IT^ . There are no interference terms (these would appear 
from the square of the sum of matrix elements): interference only contributes at the 
next order in the 1/A^ expansion. 



6 Weak coupling 

In the limit of weak coupling, the integral equation simplifies considerably and it is 
possible to obtain an approximate analytical solution. At weak coupling, we make 
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the replacements 



Dr/a{P)^-^, /^^(^)^-(^) 2Imn^(P), (124) 



so that the kernel becomes 



A 



A(/2,P)__6( — ) ImU^iR-P), (125) 



and the width reads 

2 



Tp = ^ ^ p{R) ImUniR - P) [n(r°) - n(r° - c^p)] . (126) 

We will use the integral equation in the form 

U^T^X{P) =P' + l [ Hr' - ^p) - nir')] ^^P2(p ■ r)x(r)p(i?)A(i2, P). (127) 
Jr ^ 

Both on the left and the right-hand side of this equation two integrals (over r and 
k = |r— p|) remain to be done. The imaginary part of the single bubble Imll/j is known 
analytically, see Eq. For this reason the weakly coupled case is considerably 

easier than the full problem. 

For the massless weakly coupled theory, we will now show that it is possible to find 
an approximate analytical solution in the limit p —>■ oo. Consider first the thermal 
width on the LHS. Evaluating the two remaining integrals in that limit one finds (see 
also Appendix G in Ref. [22] for the case = 1) 



lim T„ = 

P-oo P 23047r pN 



1 \^T^ r 24C(3)T 



■P/T\ 



TT^ p 



(128) 



Neglecting for a moment the ladder contribution on the RHS of the integral equation, 
it is easy to solve for x{p) ^iid we find 

lim xip) = 23047r— ^ [1 + O {T/p)] . (no ladders). (129) 

To include the ladders we now use the same momentum dependence for x but with 
an arbitrary prefactor 

X{p) = Kp\ (130) 
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Inserting this ansatz on the RHS of the integral and performing the integrals while 
consistently dropping terms that are suppressed with respect to the leading be- 
havior, we find the following result for the complete integral equation: 

Kp^ [1 + 0{T/p)] = / + [1 + 0{T/p)] . (131) 



23047rA^ ' ' ' " " ' 34567riV 
The solution for ultrahard p is therefore 



hm x{p) = mi2Ti-^. (132) 



A comparison between the asymptotic solution without ladders, Eq. ()129|) . and with 
ladders, Eq. ()132|) . reveals that the effect of summing ladders is simply algebraic: the 
full vertex V{p) in Eq. (j85|l equals 3. 

Using this asymptotic form of x{p) expression ()8(ij) for the viscosity yields 

r,^^'J^!^.9U,^,^, (133) 

where the subscript indicates that this solution has been obtained for ultrahard mo- 
mentum p. 

We can compare this result with the numerical results obtained for the shear 
viscosity in Refs. [23 ^] for the weakly coupled massless = 1 theory. In order to 
do this we must find the full dependence at weak coupling, not just the leading order 
result cx A^^. For this we use again the 2PI effective action, but now employing the 
loop expansion to three loops. The important formulas are summarized in Appendix 
A. The result is that the full A^ dependence of the shear viscosity is quite simple, and 
we find 

^ 27648C(5) N-^ 11^^0419^^ (134) 
TT A^ + 2A2 Ar + 2A2- ^'""^^ 

The numerical constant is extremely close to the ones obtained numerically in Ref. 
(3040) and Ref. ^ (3033.5). 



7 Variational solution 

In order to obtain the shear viscosity for general values of the coupling constant and 
mass parameter, we solve the problem of extremizing the functional Q in Eq. (|99|) 
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variationally. Following Arnold, Moore and Yaffe OH 1^. we expand x{p) ^ 
finite set of suitably chosen basis functions 

x(p) = iV5^a„</.(™)(p), (135) 

m=l 

where we factored out an explicit factor of N, so that the integrals below are A^- 
independent. Using this Ansatz in the functional Q yields 



(136) 



with 



/■oo 

5,„= / dp5(p)0(™)(p), 

J^mn=N dp ^(p)0(™nP)0^"Hp), (137) 

Jo 

/"OO /"OO 

n^n = N dp drn{p,r)(l)^'^\p)(l)^^\r). 



Note that iS^ is a 1-dimensional integral, J?-m„ and Timn for the line diagram are 
3-dimensional integrals, and Timn for the box diagram is a 4-dimensional integral. 
Extremizing the functional with respect to the variational parameters gives a 
simple linear algebra problem. The shear viscosity is given by 

V = -TTT^^S^am, = ^{J^ - n);^„Sn. (138) 



307r 



m 



Given the asymptotic solution xip) ~ (p/T)"^, our choice for the trial functions is 

^P'-^'iP) = (1 +''//?)li E^^^ ^139) 

For fixed A^ar these basis functions are a linear combination of the functions used 
in Refs. |191 1^ 1^ . For given basis functions the integrals in Eq. fjl37|) can now be 
performed numerically. We found that the presence of the mass M in the integration 
boundaries for p and r in Tibox (see Sec. 5.2) reduces the effort required for the nu- 
merical integration compared to the massless case. Using straightforward quadrature. 
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Figure 16: Shear viscosity rj/rjao as a function of m^/T, for four values of A(/i = T). 
Weak denotes the result in the limit A ^ 0, i.e. keeping only the single bubble in the 
rung. 



the errors due to the numerical integration are smaller than the width of the lines in 
Figs, [ini CZI Also, the results shown in these figures are obtained using three basis 
functions. Again, the effect of using this truncated basis set is smaller than the width 
of the lines. 

In Fig. ^1 we present the shear viscosity as a function of the renormalized mass 
in vacuum. To remove the trivial N'^T^/ dependence, we scaled the result with the 
approximate analytical solution rjoo (see Eq. ()133p ). The line labeled 'weak' shows 
the result in the weak coupling limit A ^ 0, i.e. preserving only the single bubble in 
the rung. The other lines represent the solution to the full problem for various values 
of A(T). We have chosen /i = T to present our results, but the shear viscosity is RG 
invariant. We observe that the shear viscosity has a characteristic dependence on the 
mass, but a negligible dependence on the coupling constant, after the dominant 1/A^ 
behavior has been scaled out. 

In order to analyze this in more detail, we show in Fig. El the shear viscosity as a 
function of the coupling constant A(T) for various values of the renormalized mass at 
zero temperature. Recall that because of the scaling with rjoo, we concentrate here on 
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Figure 17: Shear viscosity rj/rjao as a function of \{^ = T), for 4 different values of 
rriR/T. 

the subdominant A dependence. We find an appreciable dependence on the coupling 
constant only for vanishing mass. Comparison with the thermal width in Fig. 
shows that for small m/j, the shear viscosity and the inverse width have a similar 
dependence on the coupling constant, as expected. 

8 Conclusions 

We have presented a diagrammatic calculation of the shear viscosity in the 0{N) 
model at first nontrivial order in the large N limit. The 1/N expansion of the 2PI 
effective action at next-to-leading order leads in a straightforward manner to the dia- 
grams contributing to the shear viscosity at leading order and provides automatically 
the integral equation required to sum them. 

In the weakly coupled massless theory, the integral equation could be analyzed 
analytically in the limit of ultrahard momentum. Using this result we found for the 
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shear viscosity in the 0{N) model 



9216C(5) 3N^ T3 
n iV + 2 A2 



3041.9 



iV + 2 A2 ■ 



3iV3 7-3 



(140) 



This result is extremely close to the numerical values determined earlier for the single 
component theory [N = 1). 

For the general case we computed the shear viscosity numerically through a vari- 
ational approach. The results are presented in Figs. El and El Factoring out the 
basic l/A^ dependence, the results show that the remaining dependence on the cou- 
pling constant is very weak, while the effect on the mass parameter is larger. For 
the allowed range of parameters, we conclude that the shear viscosity is close to the 
result obtained in the weak-coupling analysis. 

From a more general point of view, we found that the availability of an effective ac- 
tion that sums the appropriate diagrams is extremely useful. While the computation 
of transport coefficients is still quite involved, it is streamlined by the organization 
inherent in the 2PI effective action. From the wider perspective of nonequilibrium 
quantum fields, it is satisfying that 2P1 effective action techniques can be applied 
both far from equilibrium, relying mostly on numerical tools, as well as (very) close 
to equilibrium, as we have demonstrated here. 

As we have argued recently the first nontrivial truncation of the 2PI effective 
action also sums the relevant diagrams to obtain the shear viscosity and the electrical 
conductivity in gauge theories to leading-logarithmic order in the weak coupling limit 
or to leading order in a 1/Nf expansion, where A'^^ is the number of fermion fields. A 
diagrammatic analysis of transport coefficients in gauge theories beyond the leading- 
log approximation at weak coupling is not yet available. A convenient starting point 
for such an analysis may be based on a more general nPI effective action approach j39j . 
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A Three-loop expansion 

In order to find the full A^ dependence of the shear viscosity in the 0{N) model in 
the weak coupling limit and not just its leading behavior at large A^, we summarize 
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here the results for the three-loop expansion of the 2PI effective action. We write 



XN + 2 
~8 3 
iX^ N + 2 
48 3N 



G^(x,x) 



G\x,y). 



(141) 
(142) 



xy 



Since mass and coupling constant renormalization do not enter here, we denote the 
coupling constant with A. We also immediately specialized to Gab = ^abG. The 
corresponding self energies are 



si? (^,2/) 



AiV + 2 



and the kernel reads 



^Scd(^' y') = l^abScd + SaJbd + SadSbc] 

x6cix - y)5c{x - y')5c{x' - y), 

^ab cdi^y 2/; X', y') = -T^^ [4:SabScd +{N + 6){6acSbd + ^adSbc)] 



2 3A^ 
~Q 3iV2 



A 



SabG{x,x)Scix - y), 
6abG^{x,y), 



(143) 
(144) 



18iV2 

xG'^{x,y)6cix - x')6ciy - y'). 



In particular we find that 

[x,y;x ,y 



-X' 



,N + 2 
3N 



G'ix,y)6c{x-x')6c{y-y'). 



(145) 
(146) 

(147) 



Using these expressions it is straightforward to find the dependence of the shear 
viscosity in the 0{N) model for arbitrary A^. Starting from the large result in 
the weak coupling limit, the thermal width Fp and the kernel A are related to the 
thermal width and the kernel presented in Sec. El as 



A{R,P) 



arbitrary A'' 



arbitrary A'' 



N + 2 

N " 
N + 2 



large TV 



N 



A{R,P) 



large 



(148) 
(149) 
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The dependence that subsequently appears in the integral equation can be absorbed 
in xip) as 

N 

arbitrary Af iV + Z 



large A'' 



SO the final effect in the viscosity is quite simple: 

N 



7] 



(151) 

large A'' 



arbitrary A" N + 2 

In the weakly coupled massless limit we arrive therefore at the analytical result, 



9216C(5) T3 



for the shear viscosity in the 0{N) model. 
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